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Abstract 

This work investigates some global questions about cosmological spacetimes with 
two dimensional spherical, plane and hyperbolic symmetry containing matter. The 
result is, that these spacetimes admit a global foliation by prescribed mean curvature 
surfaces, which extends at least towards a crushing singularity. The time function of 
the foliation is geometrically defined and unique up to the choice of an initial Cauchy 
surface. 

This work generalizes a similar analysis on constant mean curvature foliations and 
avoids the topological obstructions arising from the existence problem. 
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1 Introduction 



The flavour of General Relativity viewed as an initial value problem for the field equations 
is based on the geometrical nature of the equations, which implies the diffeomorphism 
invariance and the absence of a metric background structure. This has consequences for 
spacelike foliations, in that the time function of a foliation is not canonical, but arbitrary 
unless it is tied to some geometrical quantity. The latter would turn the analysis of the 
global structure of spacetimes into an investigation of the asymptotic behaviour of such a 
foliation. 

In spacetimes with certain symmetries some global foliations, defined by time coordinates 



tied to the symmetries are known, e.g. |Re], | An 



A geometrically defined foliation which does not depend on the symmetries are the well 
known constant mean curvature (CMC) foliations, where the time coordinate is given by 
the mean curvature of the leaves, which varies continuously from leaf to leaf, see |R1] for 



a survey of this topic and [RA\, [R6] for foliations of spacetimes with symmetries. 

Unfortunately the CMC constructions suffer from the existence problem, which is unsolved 
in general. To overcome this difficulty, a foliation with leaves of prescribed mean curvature 
(PMC) has been constructed at least locally in time in [H2| ] for cosmological spacetimes. 
The prescription is given implicitly, letting the mean curvature vary continuously along the 
normal vector field of the foliation relative to a given Cauchy surface. The time function 
of the foliation is geometrically defined and turns out to be intrinsic, in that coupled to 
Einstein's field equations adapted to the leaves, one obtaines Cauchy data for spacetimes 
foliated by PMC leaves. 

The aim of the present paper is to globalize this result for certain spacetimes. Motivated 



by the method in |R3|, where satisfying results about cosmic censorship have been ob- 
tained for spatially homogeneous models, I consider here cosmological spacetimes with 
two dimensional spherical, plane and hyperbolic symmetry. This choice is the first step of 
successively lowering the degree of symmetry to obtain more and more general results. I 
focus on the following questions: 

How large is the maximal interval for the time function and does this maximal foliation 
cover the whole spacetime? 



One guideline taken from [US] is, that the foliation may be extendible as long as the mean 



curvature of the leaves remains finite. We will see, that this principle stays true with some 
modifications. 

Another important aspect for the construction of global foliations is the choice of an ap- 
propriate matter model. The key requirements for the matter fields are their regularity 
in a regular geometric background as well as some energy conditions. These requirements 
are not trivially satisfied, since the first one rules out matter such as the perfect fluid, 
which is known to develop singularities (shocks) in a regular geometric background. For 



dust it has been shown in [R5\, [IR|, that there is no way to construct a global CMC foli- 
ation, which covers the whole dust-fllled cosmological spacetime. Thess counterexamples 
emphasize dramatically the importance of choosing appropriate matter models. Further- 
more demanding energy conditions seems natural and obvious, but we will see, that the 



non-negative pressures condition plays a special role in the improvements in section |3^ 



Section |^ fixes notation and states some basic definitions. We attack the main questions 



in section y, following closely the treatment in [R4]. In the final section the results will be 
discussed. 
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2 Basic definitions and formulas 



2.1 Spacetimes and foliations 

A spacetime is a pair (M,g), where M denotes a four dimensional smooth and orientable 
Lorentz manifold with metric g and signature ( — h ++). The metric induces the Levi- 
Civita connection ''V on M. If {Xq} denotes a local basis of vector fields on M, we define 
the connection coefficients "C^^ relative to this basis by '^Vxo.^p = *'C^pX^ and we get 
the Christoffel symbols T and Ricci rotation coefficients "7 by specializing to a coordinate 
basis or an orthonormal frame, respectively. The sign convention for the curvature is fixed 
by the definition ^R{X,Y)Z := "Vx^'VyZ - ^Vy^Vx^ - ^V^x,Y]Z, where X,Y Z are 
vector fields. The curvature tensor is then defined as '^R(W, Z, X,Y) := g(W,'^R{X,Y)Z) 
with Ricci tensor '^Ra/3 = "^^aiifs ^^'^ scalar curvature = '^R'll, written in abstract index 
notation of the Ricci calculus. 

Then the Einstein tensor reads GajS = ^Rap — ^^Rda/B and the field equations are Gap = 
SirTafs or equivalently "'Ra/s = Svr [Tap — ^{tr T) gap) , where Tap denotes the energy mo- 
mentum tensor of the matter fields. The matter quantities are the energy density p := 
T^i^n^n'^, the momentum density j/3 := —T^yn^h^^ and the stress tensor Sap '■= T^j^^hah'^ 
with respect to an observer, represented by a unit timelike vector n, where hap '■= 
dap + n-anf^ denotes the orthogonal projector on {n}"*- in covariant notation. 
Einstein's field equations in vacuum (Tap = 0) have a well posed Cauchy problem in har- 
monic coordinates, thus one obtains spacetimes as solutions of Einstein's field equations 
with matter, whenever the equations describing the matter fields and the energy momen- 
tum tensor couple to the field equations in harmonic coordinates, such that the Cauchy 
problem remains well posed. We will see examples in subsection |2.2| (compare [|W|, pR| ] 
for an introduction/analysis of the Cauchy problem for Einstein's equations). 

In this work we confine ourselves to cosmological solutions (M, g) of Einstein's field equa- 



tions. Due to |Ba] these are globally hyperbolic and spatially compact spacetimes, where 
the Ricci tensor contracted twice with any timelike vector is non-negative (timelike conver- 
gence condition). This last condition can be reexpressed in terms of the matter variables 
as p -|- tr 5" > for any observer, which is the strong energy condition. 

Now let us pay attention to an additional structure. A foliation {St}^ t € I (I interval 
containing zero) of (a part of) {M,g) by spacelike hypersurfaces induces on each leaf the 
unit normal vector field n, the metric hap = gap + ^anp, which also serves as orthogonal 
projection and the second fundamental form kap ■= —haK^'^V ^riy (the definition of kap 
fixes the sign conventions used in this work) . The second fundamental form is a symmetric 
tensor, intrinsic to the leaves of the foliation, and can also be written as the Lie derivative 
of the 3-metric h with respect to the normal vector field, kap = —^Cnhap- The 3-metric 
determines further geometrical objects on the leaves, such as the Levi-Civita connection 
V, the Christoffel symbols F, the Ricci rotation coefficients 7 and the curvature tensor 
R{-). Tensors intrinsic to the leaves of the foliation will carry Latin indices in the abstract 
index notation. 

The parameter t of the foliation has timelike gradient and thus can be regarded as 
(coordinate-) time. Given local coordinates (x^) on 5o, we can Lie-transport them to 
neighbouring leaves along an arbitrary family of transversal curves, parametrized by x. 
We will express equations containing coordinate components with respect to the adapted 
coordinates (t, x). 

The lapse function N and the shift vector 1/ _L n on the leaves are defined by dt = Nn + v, 
thus A'^ = —g{dt,n) and v = dt — Nn. Then we have 1 = dt{dt) = N dt{n). Further, dt 
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is (co-)orthogonal to the leaves and if we denote the conormal of the leaves by a we see 
that dt = —N~^a or o" = —Ndt, thus measures the length of dt. 
The most common example involving the lapse function is the event horizon in Schwarz- 
schild spacetime, where the coordinate time explodes along the worldline of infalling ob- 
servers n, thus explodes or A'^ — > 0. In this work we will be faced with a com- 
plementary scenario, where we have to ensure, that the lapse function does not explode, 
corresponding to the phenomen of recollapse, where coordinate time freezes, as dt — > 0. 
If this occurs, one could try to reparametrize the foliation by setting t = f{t) for some 
monotone function / and one gets dt = f'{t) dt, 

and the same relation holds for shift. 

The 3 -|- 1 split of the spacetime geometry by means of lapse and shift ends up in the 3 + 1 
form of the field equations. The constraint equations are 

(la) R + — |A;p = IGvrp (Hamiltonian constraint) 
(lb) y-'kij — ViH = Sirji (momentum constraint) , 

with |A;p = kafjk'^^ and H = trk denotes the mean curvature of the leaves. The ADM 
equations read 

(2a) dthij = -2Nkij + Vii^j + Vjfj 

(2b) dthj = -ViWjN + iV (^Rij + Hkij - 2klkrj - 87r( Sij + \{p-ii S)hij 
+i''''Vy.kij + kj-jViV^ + kirVjV^ 

Taking the trace of the second equation and eliminating the scalar curvature R by the 
Hamiltonian constraint we obtain the lapse equation 

(3) AN + N(^\k\^ + A7:{p + tvS)'^ ={dt-u)H , 

which serves as a constraint of the foliation. Note, that in cosmological spacetimes, the 
term in brackets is always non-negative. If it turns out to be non-vanishing, then the 
left-hand side of (^) can be shown to be an isomorphic mapping of N, considered as an 
element of some Sobolev space H'^ into H^~^. This observation motivates the following 

2.1. Definition. 

A Prescribed Mean Curvature (PMC) foliation is defined to be a foliation {St} 
satisfying (|3|), with 

(4) {dt - v)H = Nn{H) := 1 , 

thus the mean curvature of the leaves is forced to vary uniformly along the normals 
of the leaves. 



In [tI2] I proved the following local in time result: 
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2.2. Theorem. 

Let {M,g) be a smooth, globally hyperbolic spacetime, obeying the strong energy 
condition, with compact Cauchy surface S and 

(5) A = +47r(p + tr5) 

does not vanish identically on S. 

Then there exists aT > and a unique smooth PMC foliation {St}, t € [— T, T] in 
{M,g), with S = S'o. 



Note, that the setting here is quite general, no symmetry assumptions have to be made 
and essentially the strong energy condition turns out to be sufficient for the local in time 
existence of a unique PMC foliation up to the choice of an initial Cauchy surface. 

The aim of the present paper is to globalize the result. Here are two problems involved: 
How large is the interval of values taken by the time coordinate and does the global 
foliation then cover the whole spacetime? To answer these questions in general there seem 
to be no techniques available up to now. One strategy to obtain global results is, to study 
first spacetimes with some spatial symmetry, taking advantage of the simplifications of 
the equations. Then the hope is, that the techniques developed in these cases give insight 
into the nature of more general classes of spacetimes, by successively lowering the degree 
of symmetry. Here we will focus on spacetimes with two dimensional spacelike orbits of 
symmetry and three (local) Killing vector fields. In the second part of this paper we 
will consider spacetimes with two commuting (local) Killing vector fields. These cases 
are indeed the first steps of this program, since the initial analysis for locally spatially 
homogeneous spacetimes has already been done and lead to very strong results (see 
for the exact analysis or [HI] for an overview about the results in the present context). 



2.2 Matter models 

Before getting deeper into the analysis just motivated I introduce some matter models and 
their coupling to the Einstein equations, with special emphasis on the Cauchy problem. 
Since energy conditions will play an important role in the estimates we perform later on, 
I assemble the relevant ones for the present work first. 

• The dominant energy condition. Its statement is, that for all orthonormal frames 
{^o} with ''eo timelike 

TCeo,'eo) > \TCe^,'ep)\ . 

Another formulation of this condition is, that for any observer, the local energy 
density p is non-negative and the momentum density j is non-spacelike, thus 'matter 
cannot travel faster than light', a statement, that can be proved rigorously (see for 
example section 4.3 of pIE[ ), leading to the result, that if the energy momentum 
tensor obeys the dominant energy condition and vanishes on a set S, then it vanishes 
on the whole Cauchy development D{S) of the set, thus D{S) is a vacuum spacetime. 

• The strong energy condition states that for all timelike vectors v the inequality 

Tiv,v) > l{tTT)g{v,v) 
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holds, or equivalently for any observer p + tvS > 0, thus the stresses do not become 
too negative. Another formulation is '^R{v, v) > 0, also known as the timelike conver- 
gence condition, which contributes to the expansion of timelike geodesic congruences 
a negative term, thus shifting the balance towards contraction to a final singularity 
of the congruence. Hawking's famous singularity theorems guarantees then geodesic 
incompleteness in the past provided further the existence of a Cauchy surface with 
uniform negative mean curvature. 

We already used the strong energy condition as an integral part in the definition of 
cosmological spacetimes and its meaning for foliations. 

• The non-negative pressures condition demands the stress tensor S to be positive 
definite. This condition ensures in some sense, that the pressures contribute more 
to attraction than to repulsion, leading to a finite lifetime of the spacetime under 
certain circumstances. This somewhat unexpected behaviour is a true relativistic 
effect, which does not have a Newtonian counterpart. 

2.2.1 The Maxwell equations 

The Maxwell field is described by a two form F, subject to the Maxwell equations 

(6a) dF = ^ '^[.F^u] = 

(6b) d*F = 47r*J ^ " V^F"/^ = 47r J" , 

where J denotes the electromagnetic charge current density. Alternatively the covariant 

derivatives can be replaced by ordinary derivatives. 

It is well known, that the Maxwell equations admit (locally) a reformulation in terms of the 
vector potential A, a one form with dA = F. Then the remaining inhomogeneous equation 
reads "V^^VqA^ — ^V^^V^Aq, = 4TTJa- Fixing the gauge invariance of the equation by 
the Lorentz gauge, we obtain the system 

(7a) * = (Lorentz-gauge) 

(7b) nAa = 4TrJa - RatJiAl' , 

with □ = — " V'^'^Vju and the curvature term arises as a consequence of some commutations 
of derivative operators. In this formulation the second equation is a wave equation, hence 
can be brought into first order symmetric hyperbolic form, if the source term is appropriate. 
Thus we get local existence and uniqueness for this equation in a given spacetime. One 
can further show, that then the Lorentz gauge propagates and we get indeed a unique local 
solution of the Maxwell equations. Spatially global solution can be obtained by the usual 
patching argument, arising from localizing the equation with respect to an appropriate 
partition of unity of the initial data on some Cauchy surface. 

Given an electromagnetic field F we can form the associated energy momentum tensor E, 
defined as 



(8) Eap = ^ (Fa^F/ - l\Ff g^p) , 
where = F'^'Fij,^. 

E is tracefree, satisfies the relation *V^-E'°'* = 
energy condition hold. 



—F°'nJ^^ and the dominant and strong 
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Since the Maxwell equations turn out to be symmetric hyperbolic at least locally, we can 
couple them to the Einstein equation to get a symmetric hyperbolic system of equations in 
harmonic coordinates. Thus we end up with a well-posed Cauchy problem for the Einstein- 
Maxwell system, as long as the electromagnetic charge current density is in appropriate 
form. 



2.2.2 The Vlasov equation 



The Vlasov equation is a model for a collisionless gas. It describes the motion of a huge 
number of structureless particles in spacetime. We need only the case, where the particles 
have unit mass, where the equation is composed of a non negative function /, defined on 
the mass shell of particles of mass one P := {v (z TM | g{v,v) = —1, future pointing}, 
representing the particle distribution, and a geodesic spray X on P. The equation reads 

(9a) X{f) = , 

with 

(9b) ^=P^5,-^r^.pV^ = ^"^e«-^7^^^^^^ , 

where denotes the components of the momentum of the particles with respect to the 
given coordinates and are the components of the momentum with respect to an or- 
thonormal frame. They are related by = (^Cq,)'^ and on the mass shell is determined 
by the components and = \/l + 6ijV^v^ . 

Inserting the definition = (%a)^9^ and doing a 3-1-1 decomposition (with the exception, 
that we do not write down the explicit expression for *7oj) we can reformulate the Vlasov 
equation as 

= dtf + {N${e,r-i^^)dJ 

- (e,(iV)^o + Ni-krMYie.r + 6,,%y + N6,kl^^-^) ^/ ■ 



The Cauchy problem for the Vlasov equation in a given spacetime is easy, since it is a 
linear, scalar equation with characteristic vector field (Y, Q) satisfying 

and since X{f) = 0, / is constant along the characteristics. 

From the particle distribution / one can construct other physically meaningful quantities 
such as the energy momentum tensor by integration over the tangent spaces. We denote 
the part of the mass shell in the fibre over x ^ M hy P^ := P DT^M . Then we define the 
energy-momentum tensor T by 

(10) Tap{x):=- fPaPpV\9\/Podp^dp'^dp^ ■ 

T is divergence free, satisfies the dominant and strong energy condition and the non- 
negative pressures condition. 

To obtain the matter quantities p, j and 5 one has to calculate the components of the 
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energy momentum tensor with respect to an orthonormal frame. These components have 
the following representation: 



Now let us consider the Einstein- Vlasov system. The energy momentum tensor automati- 
cally satisfies divT = as mentioned above, thus there arise no additional equations from 
the Bianchi identities. But unfortunately the coupled system of equations is not sym- 
metric hyperbolic in any sense, due to the fact, that it is a system of integro-differential 
equations. Nevertheless, the local existence proof applies similar techniques as in the case 
of quasilinear symmetric hyperbolic systems where the peculiarity in the construction for 
the Einstein- Vlasov system consists in bounding the support of / in the tangent space: 
There is no a priori bound on the velocities of the particles, and no localization argument 
available as for the spacetime coordinates. In order to estimate the matter quantities 
appearing in the coupled system, one has to control the maximal velocity uniformly dur- 



ing the construction. This has been done for example in |R2|, establishing a well-posed 
Cauchy problem for the Einstein- Vlasov system. 

With this result at hand, it is easy to extend this result to the Einstein- Vlasov-Maxwell 
system, since the Einstein-Maxwell equations are symmetric hyperbolic and when coupling 
the Vlasov equation to them nothing new appears. 



2.3 Further conventions 



For the convenience of the reader I cite the well known Gronwall estimate, which is central 
to the analysis of partial differential equations and will be used here in this work. I adopt 
the convention of denoting any generic constant by C. 

2.3. Proposition (Gronwall's inequality). 

Let I cM. be an interval, to G I and a,(3,u G C(/,M+), with 



u{t) < a{t) + 



I3{s)u{s) ds 



to 



for all t & I. 
Then 



a 



(s)/3(s) e\^s^^^)'i^\ ds 



to 



u{t) < a{t) + 
holds for allt^I. 

The proof of Gronwall's inequality in this particular form can be found in . 
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3 Spacetimes with two dimensional spherical, plane and hy- 
perbolic symmetry 



3.1 The geometry of surface symmetric spacetimes 

Let (M, g) be a smooth, globally hyperbolic spacetime which is topologically of the form 
X X F, with F a compact, orientable surface. The submanifolds {r} x x F are 
assumed to be Cauchy surfaces of M. The universal covering F of F induces a spacetime 
(M, ^) by M = M X 5^ X F and g = p*g, p : M — > M the canonical projection. Moreover, 
there is a group G of isometries acting on (M,g). 
Then (M, g) 

• is called spherically symm,etric, if F = and G = SO{3) acts isometrically and 
without fixed points on x S*^ 

• is plane symmetric, if F = and G = E2 (Euclidian group) acts isometrically on 

F = M? 

• has hyperbolic symmetry, if F has genus greater than one and the connected com- 
ponent of the symmetry group G of the hyperbolic plane, H^, acts isometrically on 
F = (thus F = H^/T, with T a discrete group of isometries of H^) 

and the matter quantities remain invariant under the isometries. 

To collect these cases, each such spacetime is called surface symmetric, the diffeomorphic 
images of F in the product decomposition of M surfaces of symmetry and each surface in 
M diffeomorphic to x F will be called symmetric. 

In expressions involving indices, lower case Greek indices range from to 3, lower case 
Latin indices (preferably taken from the middle of the alphabet) range from 1 to 3 and 
upper case Latin indices (from the beginning of the alphabet) take the values 2 or 3. 

The isometric action forces the curvature of the surfaces of symmetry up to rescaling to be 
e = 1,0, —1 in the spherical, plane and hyperbolic case, respectively. Therefore they can 
be coordinatized by the well known angles tp) which cast the metric g of the surfaces 
of symmetry (considered for a moment as abstract manifolds) into the form 



' sin I? 6 = 1 

1 , e = 

sinh'f? e = —1 



(11) g = di^'^ + el dif"^ , := < 
Define the area radius function r on a surface of symmetry F (embedded in M) to be 



(12) r = ^^Yol{F) , 

then r is independent of (i!), if) and the metric of F reads 

(13) g = r^~g . 

With respect to any symmetric Cauchy surface S we have the timelike unit normal vector 
n of in M. Regarding n as a normal vector to F in M, we can define a second unit 
normal m of F by the conditions, that m is tangent to S and the system (n, 771,62,63) 
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is positively oriented, where we set ba '■= r ^ca, ^2 '■= 63 := d^,. We denote the 
associated second fundamental forms of (F, g) in (M, g) by n and A, with 

(14a) \ab = -\m{gAB) = -\'rn{r'^) qab = \{tTX)gAB , trA = -^m{r) 
(14b) KAB = -^n{gAB) = -^n(r^) gAB = ^(tr k) gAB , tr k = -^n{r) 

Consider now a Gaussian coordinate neighbourhood {x',i!),ip) of a surface of symmetry 
F, covering (a part of) a symmetric Cauchy surface S. The metric h of S then takes 
the form h = dx'"^ + \/\}i\g. The projection of geodesies starting in (M, g) orthogonal at 
F remain orthogonal to all surfaces of symmetry. Following them until their projection 
meets F again, the symmetry allows only two possibilities: The point of return is the 
same as the starting point or an antipodal point, in which case we force the geodesic to 
turn a second time around the circle. Let L denote the length of the geodesic. Setting 

a = 2tt (^jl^\h{z)\~^/^ dz^ , we define a new coordinate x{x') by a Jq \h{z)\~^/^ dz. In 

the coordinates (x*) = {x,'d,ip) the metric has the representation 

(15) h = {dx"^ + a'^g) = A'^dx^ + g , 

with A{x) = a~^\h{x)\^^^ defined on S^. Comparing this with equation (|l3|) shows r = Aa. 
The corresponding Laplacian A acting on a function ip on S can now be calculated to 

= -A-2(V;" + A-^^V) + , 

where A denotes the Laplacian of g, and the prime differentiation with respect to x, a 
convention, which we adopt for the rest of this work. 

Furthermore, the symmetry and the given coordinate representations permits the second 
fundamental form k to have the form 

(16) k = A'^Kdx'^ + l{trK)g , 

where the coefficents are functions on S^. Taking the trace yields the mean curvature 
H = ilk and we get the relation 

(17) H-K = tTK . 



So far we know the intrinsic and extrinsic geometry for a symmetric Cauchy surface S in 
M. Let us turn now to the 3+1-geometry. Theorem 2.2 states the conditions to guarantee 
local in time existence for a PMC foliation {St} of a neighbourhood of S. We need, that 
A defined by equation (^) is non- negative and does not vanish identically on E. Again n is 
the unit normal on S and T the energy-momentum tensor of the spacetime. Therefore we 
get a sufficient condition by the assumption that M satisfies the strong energy condition 
(hence the matter term in (^) is non-negative) and that there exists at least one point 
in E, with A > 0. The strictness of the inequality is not a restriction, since one always 
can perform a slight deformation of S, such that the second fundamental form is not 
identically zero, which does the job. The same reasoning works, of course, if M satisfies 
the dominant energy condition and the non-negative-pressures condition (then both parts 
of the matter term in (|5|) are non- negative). Thus, in surface symmetric spacetimes, some 
energy conditions are sufficient for the existence of a local in time PMC foliation. 
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Having constructed a local in time PMC foliation, one can ask, if the leaves of the foliation 
are symmetric, when S has been chosen to be symmetric (Note, that due to symmetry 
A = X{x) is a function on only). As described in |H2(| the PMC foliation is given as 
a limit {w,N) of functions {w^,N^) on S in some Sobolev space. Here describes a 
family of spacelike hypersurfaces in M and A^-' converges towards the lapse function of w. 
(w^ , ) are defined as solutions of the sequence of symmetric-hyperbolic elliptic systems 

with representing E and A^*^ = 1, hence respecting the symmetries. The underlying 
metric structure of the system is the sequence of first fundamental forms /ij-i of the 
surfaces 'w^~^. If all the {w^,N^) respect the symmetries, then the PMC foliation also, 
hence it suffices to show, that {w^ ,N^) respects the symmetries, if w^~^ , N^~^ does so. 

This is clear for , because the symmetric hyperbolic equation can be localized, and the 
pullback to M is invariant under the action of the isometry group. 

Suppressing the index j — 1 from some quantities determined by the metric /ij-i, the 
Laplacian reads 

A{w^~^)N^ = -A-^ {{N^y + A~^A'{N^y) + A{w^-^)N^ . 

The elliptic equation then becomes 

((ATi)" + A-^A'iN^Y) + X{w^-\Dw^-^)N^ + A{w^-^)N^ = 1 . 

If = \/\^ = r^e^ denotes the volume form of the surfaces of symmetry F in M, one 
finds Vol(F) = /[o,,]^[o,2.]^ = 4vrr2, as desired. Setting L := - A-^ + A'^ A' £) 
and N^{x) := /[o,^]x[o,2^] = h /[o,^]x[o,27r] ^^e^' t^^^n integration of the elliptic 

equation, /[o,^]x[o,27r] + A + A)W = /[o,^]x[o,2^] = yields 

(L + A + A)A^' = 1 , 

since on the one hand /[q 7r]x[o 27r](-^ ^)^''^^ = 47r(L + A)A-^ and on the other hand 
/[o,.]x[o,2.](^^'>^ = /[o,.]x[o,2.](^^')^ = /f vanishes as well as ANK 
Thus A^-' is a solution of the elliptic equation. Uniqueness then gives us = , thus 
A^-' = N^{x), which expresses the symmetry of A^-', as desired. 

Thus we have shown, given a symmetric Cauchy surface S in M admitting a local in time 
PMC foliation, that all leaves of the foliation are symmetric, too and coordinatizing them 
in the way described above yields: 

3.1. Proposition. 

Let (M, g) be a surface symmetric spacetime obeying the strong energy condition. 
Then there are coordinates (x^) = {t, x, i?, (/?) adapted to a local in time PMC foliation 
{St} of a neighbourhood U =]ti,t2[ x5] of T, = Sq in M , which cast the metric into 
the form 

(18) g = -N'^dt^ + A^ {{dx + vdtf + a^g) . 

All coefficients except a are functions on ]ti,t2[ x 5*^, whereas a depends only on the 
time function t. A and a are everywhere positive, N denotes the lapse function of 
the foliation and v the non-vanishing component of the shift vector, uniquely fixed 
by the condition i^{t,0) = 0. 
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Supplementary to the notation already introduced, let an overdot denote differentiation 
with respect to t, while keeping the prime as a marker for differentiation with respect to 
the coordinate x. 

Finally we find for the orthonormal frame {^e^} and its dual {V^} canonically induced by 
the 3 + 1-split: 

^eo = n = N-^ {do - udi) V° = Ndt 



A-^di ^a^ = A{udt + dx^) 



%i = m 

% = {Aa)-^d2 ^a'^ = Aadx^ 

"es = {Aaei))-^dz = Aae^) dx^ 



3.2 The field equations 

Given the PMC foliation, we can write down the field equations in the 3+1-representation. 
The symmetries suggest to represent the matter quantities completely with respect to an 
orthonormal frame, thus we define j := —T{n,m) = A~^ji and Sij = T(%j,%j). 
Then one calculates the constraint equations ([l|) to be 

(19) (Ai/2)" = 1^5/2 ^H' -^{H-Kf-K^- levrp) + lA'/^a'^ 

(20) K' = -3A-^A'K + A-^A'H + H' + S-nAj 

The foliation is fixed by the lapse equation (^) and the PMC condition (§) 

(21) N" = -A-^A'N' + A^n[\{H - Kf + + Att{p + tr 5)) - A^ 

(22) H = l + vH' 

and the evolution equations @i read 

(23) A = -NAK + Au' + A'v 

(24) d = -\Na{H -2,K) - au' 

(25) k = vK' - A-'^{N" - A-^A'N') 

+N (-4A-^/2(^i/2y/ ^ (^-2^/)2 8vr(^-25n + \{p - tr 5))) 

Integrating the equation for a over the circle yields d = /io, with /i = — ^ J^i N{H — 3i^), 
since J^i u' vanishes. Inserting this back gives an equation for shift: 

(26) v' = -\N{H -2,K) + \ f N{H-3K) 

Differentiation of the equation for H with respect to x yields an equation for H': 

(27) (dt - vd^)H' = v'H' 

In summary we have equations for space and time derivatives of the fundamental forms h 
and k. Moreover there are equations for the spacelike derivatives of lapse and shift, but 
unfortunately there is no information about their time derivatives. 
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So far we started with a given surface symmetric spacetime, admitting a local in time PMC 
foliation, to which we adapted the field equations. This raises the opposite question: Given 
the equations ([l9|)-(|27|) and appropriate data, does there exist a solution, and how unique 
is it? To answer this question we first need to state more precisely the term 'appropriate 
data': 

3.2. Definition. 

A symmetric initial data set is a smooth collection (S, h, k) consisting of a 3-manifold 
S diffeomorphic to x F with metric h and a symmetric tensor field A; on S, where 
S admits coordinates, such that h and k could be written in the form shown in 
section 3A_. 

If there are matter fields present, then it is assumed, that there is also smooth 
symmetric matter data and equations, leading to a well posed Cauchy problem of 
the reduced field equations in harmonic coordinates. 



The smoothness of the quantities appearing in the definition is required, since the trans- 
formation to harmonic coordinates involves derivatives. 

3.3. Proposition. 

Let (S, h, k) be a symmetric initial data set, with matter obeying the strong energy 
condition and ( compare equation (^) for definition ) X > somewhere on S . Further, 
let tQ denote an arbitrary real number. 

Then there exists a (5 > and a PMC foliated surface symmetric spacetime (M, g) 
diffeomorphic to ]to — 5,tQ + 5[xTi with an embedding l : S — > M, satisfying 
= StQ and L^h, t^k are the first and second fundamental form of St^ in {M,g). 
(M, g) obeys the strong energy condition and g can be written in the form described 
in proposition \3.i\. This construction is unique up to the choice of tQ and 6. 



Proof. 

On the induced manifold E diffeomorphic to xF the induced data is invariant under the 
group action, hence the Cauchy developments also and we get a surface symmetric Cauchy 
development of the data, admitting a symmetric local in time PMC foliation near S on 



some time interval ]to — d, tQ + 6[ and allowing a set of coordinates stated in proposition 3.1 . 
The uniqueness property stated above follows from the geometric uniqueness of solutions 
of Einstein's equations associated with the uniqueness of the PMC foliation, once the 
remaining degree of freedom has been fixed by the requirement = St^- Note, that 
lapse is fixed by the ellipticity of the respective PMC equation and shift is fixed by equation 
@ and the condition i^{t, 0) = 0. ■ 



In surface symmetric spacetimes, there is another way to express the constraint equations 
(|l9| ) and ( |20| ) in terms of 'optical scalars', see []GM | for an enlightening presentation. 
The geodesic null congruences determined hy k± = m it n give rise to the null expansions 



(28) 



-(tr A ± tr k) = -(m(r) ± n(r)) = -k±{r) = 2A"M' ^ [H - K) 



where (14), with (^) has been used, together with the relation r = Aa for the area radius 



defined in (12). The formula illustrates the definition of r as a volume measure, whose 
variation along k± is descibed by the (negative of the) trace of the second fundamental 
form associated to this direction. 
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Now we can write the constraint equations for an arbitrary symmetric Cauchy surface S 
symmetrically as 

87r(/j + i) = -m(i9_) - + 'd.H + r'^e 
Snip - j) = -m{i9+) - - + r'^e 

and taking ui± := rd± as the fundamental variable we get (compare |GM| ) 

(29) m(cj±) = -87rr(p^j)T^±^^ + ^(w+W--2a;| + 4e) . 

Furthermore, the area radius serves as a warping function in the warped product M = 
B Xj. F of the two dimensional spacetime {B,g) with {F,g), where B = W x is the 
quotient M/G and g = g\B- We adopt the convention that lower case Latin letters from 
the beginning of the alphabet range from to 1, and objects intrinsic or orthogonally 
projected to {B,g) will be marked by an underbar. Einsteins equations in this framework 
can be considered as equations in {B,g) for the field r : B — M: 

M 

(30) YaYb r = ^gab - 47rr {Tab - tr Tgab) , 

where T denotes the projected energy-momentum tensor into the spacetime {B,g) and 
the mass function M is defined as 

(31) M := (e - VVVar) = {e - *V"r "V^r) = Ir (e - iri?+ri?_) , 

since i^+'f?- = "'V^r^Vor. M turns out to be the Hawking mass niH^F) (up to a factor 
— \x foi' genus(F) > 3): 

(32) ^u{F) := ^^^^ (vrx(F) - 1/8 ^^^_) 

and differentiation of ( |3TD yields the mass flux equation 

(33) Ya^ = ^T^r'^{Tab-iTTg_ab)Y^r . 



3.2.1 Expanding and recollapsing models 



Let us first consider the definition of mass in (31). In the spherically symmetric case we 
see, that gradr is spacelike as long as 2m < r holds, a condition we are familiar with in 
connection with the Schwarzschild spacetime. 

In the plane and hyperbolic case the picture is quite different, since gradr turns out to 



be timelike, as long as 2m/ r > e. In fact, the lemmas 2.3, 2.4 and 2.5 in [E14| prove for 
spacetimes with hyperbolic symmetry and in non-fiat plane symmetric spacetimes gradr 
is timelike, provided the dominant energy condition is fulfilled. Thus we cannot think 
about r as some radial, spacelike coordinate any longer and this fact will play a central 
role in our further analysis: 

We can choose the time orientation in those spacetimes to find gradr past pointing, 
throughout the whole spacetime. Then we define the time orientation on the cotangent 
bundle by metric transport from the tangent bundle, such that dr turns out to be future 
pointing. Therefore r increases with time, which means, that the area of the surfaces of 
symmetry increases with time, and the spacetime expands in this sense. Thus we will 
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call these spacetimes expanding. We get an equivalent characterisation by the relation 
between gradr and the null expansions: and have fixed and opposite signs (in 
particular i}^ > 0, 'd^ < 0), since i}± = {2/r)k±{r) = {2/r)dr{k±) is the contraction of a 
timelike with a null vector. 

Therefore we can decide, given a symmetric initial data set (S, h, k) with non-flat plane or 
hyperbolic symmetry, with matter obeying the dominant energy condition, which direction 
is expanding or contracting. This is important, since we expect some singularity towards 
the contracting direction and therefore we will pay attention to the past development 
which represents the contracting direction with respect to our conventions. 

Note, that we obtained this information without referring to the mean curvature H = lik. 
But it turns out, that the mean curvature of an arbitrary Cauchy surface S in expanding 
spacetimes is also restricted in some way: 

First remember, that H = ti k = — divn measures the convergence of the geodesic con- 
gruence, future pointing and orthogonal to S. Thus, H < everywhere on S corresponds 
to the notion of expansion and Hawking's theorem proves the existence of a singularity 
in the past. In an expanding spacetime in our sense, H is not necessarily confined to be 
negative everywhere. But we will see, that it is impossible for H to become non-negative 
everywhere on S: The explicit formula ( p8| ) for '&± shows in connection with our sign 
conventions, that H < K. If H were non-negative, then \H\ < \K\, thus — < and 
integration of the Hamiltonian constraint (|l9|) over yields a contradiction. 
With the dominant energy condition we can show more. Writing the Hamiltionian con- 
straint as i? -f- — \k\^ = 167rp we get R + > 0. Assuming H = on S, we would get 
-R > 0. But symmetric surfaces S in spacetimes with plane and hyperbolic symmetry obey 
the topological conditions (i), respectively (ii), of theorem 5.2. in |SYl , which imply, that 
S cannot have positive scalar curvature and must be flat in case of non-negative scalar 
curvature. Thus 5 must be flat and the scalar curvature vanishes. This in turn forces 
p = and A; = by the Hamiltonian constraint. In the plane symmetric case then the 
spacetime is flat by /c = 0, p = and the dominant energy condition, contradiction. In 
the hyperbolic symmetric case integration of equation ([l9|) yields a contradiction, too. 

Putting all this together we conclude without loss of generality, that surface symmet- 
ric spacetimes, obeying the dominant energy condition, which are plane symmetric and 
not flat or have hyperbolic symmetry are everywhere expanding, with dr timelike future 
pointing, i!)^ > 0, "i?- < and any symmetric Cauchy surface S is not maximal with mean 
curvature not everywhere positive on S. 

Of course, these arguments do not work in the spherically symmetric case, where grad r 
is spacelike in {2m < r} and no fixed sign of the expansions ??± can be expected, fitting 
into the the general belief in the closed universe recollapse conjecture, which precludes 
expansion of the whole spacetime. In particular we expect the existence of a maximal 
hyper surf ace in M. 

3.3 A priori estimates for the field equations 

Our aim now is to get sufficient estimates, that allow the construction of a global PMC 
foliation. The 'size' of the foliation is measured by the mean curvature H, thus we are 
looking for uniform estimates of the geometric and matter quantities in terms of H. 

So let (M, g) be a surface symmetric spacetime, satisfying the dominant and strong energy 
condition. Let us assume, that there is a Cauchy surface S in M with mean curvature 
H < 0. In particular we get a local in time PMC foliation {St}, t G t2[ with S = 5*0 by 
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theorem 2.2. If the spacetime possesses plane or hyperbolic symmetry, we choose the time 



orientation in correspondence to the conventions introduced in 3.2.1, thus H decreases with 
decreasing PMC time. If {M,g) is spherically symmetric we choose the time orientation, 
that H decreases with PMC time, too. Then in either case we expect to find a singularity 
at least in the past D'(T,) of S. 

In D~{'E,) the mean curvature is bounded from above hy H < H < and H = H only 
on S. Thus is bounded from below and we find the following estimates for the field 
equations in D~{T,) as long as H remains finite: 



At first we consider the constraint equation (|29D on a fixed leaf. At the critical points 
of co± we find together with the dominant energy condition the important inequality 



(34) |ri?±| < 2(\Hr\ + v^(i?r)2 + e 



as shown in R4|. 



For plane and hyperbolic symmetry this inequality ( p4D can be strengthened to 

|i?dz| < MH\ <C =^ \A-^A'\ < C and \K\ < C 



by the definition ( |28D of the null expansions. 

In the spherically symmetric case the argument is more complicated: First, the work 
Q of Burnett shows, that under the additional assumption of the non-negative- 
pressures condition we have r < C and < C < M. Using the upper bound for r on 
the right-hand side of (|3^ ) we get with (|3l| ) 

|n?±| <C =^ - <C 
r 

Thus we get < Cm~^ < C, hence r is bounded above and below away from zero. 
Inserting this in the estimate (|3^ leads to < Cr~^ < C and we are in the same 
situation as in the plane/hyperbolic case, hence |^~2yl'j < C and \K\ < C. 

Now consider the lapse equation (^ij) on a fixed leaf. At the point, where N attains 
its maximum N, we have 

N < (^^{H - Kf + K^ + Att{p + tr S)) < C/H^ < C , 

due to the strong energy condition. Hence we have a bound for |A^|. 



• Next, there is a bound for shift. Examination of formula (|26| ) shows, that all quan- 
tities on the right hand side are bounded, so we get < C, and therefore, using 
i/(t,0) = 0: \u{t,x\ < \i^{t,0)\ + fg^\iy'\ < C. 



By the way, the equation ( |27D for H' provides a bound for \H'\, since the coefficient 
on the right-hand side is bounded and applying Gronwall gives the desired estimate. 



With the information about u' examination of (24), a = a (^ — ^N{H — 3K) — u'^, 
shows, that the factor in brackets is already bounded and we get an inequality of 
the form |5t ln|a|| < C, which leads to a bound for \a\ and |a~^|. 



The same line of argumentation works for A. Equation ( p3D can be written as 
A — vA! = A{—NK -|- z/'), with the factor on the right-hand side bounded. So we 
get bounds for \A\ and too. 

Since we have already bounded the null expansions t?± = 2A~^A' it {H — K), one 
sees easily, that even \A'\ < C. 
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Integration of equation (|19|) over the circle yields an inequality 

< I I A^'I^H^ + \ f A'/'^a-^ . 

From this one concludes the boundedness of J^i p, and by the dominant energy con- 
dition the boundedness of Jgi\j\ and /^il'S'l. Integrating now equation (|2ll) starting 
at a point, where A^' = 0, we get 

= - J A-^A'N' + J A^n(^^{H - Kf + + A7r{p + trS)^ ~ j 
\N'\ <C+ f \A-^A'N'\ . 



The bound for |A^'| then follows from Gronwall's inequality. 

Furthermore, the bounds for A,A~^ together with the basic estimate for ■i)± and the 



boundedness of J^i p are enough to apply the proof of the lemma in [R,5], which ends 
up with |iV~^[ < C. 



Collecting all these estimates we get the 
3.4. Proposition. 

Let (M, g) be a surface symmetric spacetime, obeying the dominant and strong energy 
condition and in the spherically symmetric case the non-negative-pressures condition, 
too. Assume the existence of a symmetric Cauchy surface S with strictly negative 
mean curvature. In particular we get from proposition |3. j| a PMC time coordinate 
t, ranging in ]ti,t2[ with S = {t = 0} and H decreases with decreasing t. 
Then we have uniformly on ]ti,0] 

|yl|, |yl"^|, l^'l, |a|, |a"^|, \H\, \H'\, \K\, \N\, |iV-i|, |Ar'|, < C . 



To put this result into some framework, we establish some formalism, to have useful abbre- 
viations at hand as well as to make clear the dependence between estimates of geometric 
quantities and matter variables. 

3.5. Definition. 



r := [A,a,N, i^,H,K) 
collects the quantities describing the geometry of the foliation and 

ip,j,S) 

abbreviates the matter quantities. 



We have already estimated the quantities J-, A~^ , a^^ , N^^ , as well as A' ,N' .,u' ,H' and 
A,d,H (by inspection of the field equations). The idea is now, to bound all quantitities 
J^, ^ together with all of their derivatives uniformly on the time interval 0]. Then there 
exists a smooth extension to the closure of the interval, which serves as a new symmetric 
initial data set for the field equations in the sense described in definition 3.2. Note, that the 
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bounds for a-\ A^, A^"^ ensure, that the geometry remains regular at the boundary 
of the time interval and the C°°— bounds for lapse and shift turn out to be necessary to 
obtain C°°— bounds for the fundamental forms. 



Proposition 3.3 then sets us in the position to extend the foliation at least in the past 



direction, where H < holds: Construct the solution stated in proposition 3.3 and embed 
M into the maximal Cauchy development of the data. 

To carry out this program, we need some knowledge about the matter quantities. First, 
the matter has to obey the dominant and strong energy conditions. Second, assume that 
the regularity of the geometry guarantees the regularity of the matter in a certain way: 
For all non-negative integers m and n we have 

\drd^T\ <C =^ \drd^^\ < C 
^^^^ \drd^+^T\ <C =^ |5r^^9"^| < C , 

then the following lemmas hold. 

3.6. Lemma. 

Assume, that the matter fulfills p5|) . Then for arbitrary non-negative integers m,n 



Vk„ \dTdiT\ < C 



\drd^+^T\ < C 



holds. 



3.7. Lemma. 

Assume, that the matter fulfills (|35|) . Then for arbitrary non-negative integers m,n 



Vk„+i \drdiT\ < C 



d^^^d^Tl < C 



holds. 



Together with proposition 3.4 the lemmas accomplish the task of bounding all neces- 
sary derivatives of the geometric quantities, by first bounding all spatial derivatives d'^J- 
(proposition 3.4 and lemma |3^ and then successively all derivatives d^d^^J-., with k-\-l = n 



by alternative applications of the lemmas. Therefore, in view of proposition |3.4| the valid- 
ity of property ( |35| ) will be enough to extend the foliation. We will prove ( |35[ ) for some 
matter models in the next section. But first, of course, we have to prove the lemmas. 

Proof of lemma W. 



A: Applying on equation ( [T9|) and integrating along yields a bound for the 

difference \{dY'd'iA^I^)'{y) - {d'l'd!iA^/^)'{x)\ and using Jg^idY'd^A^/^)' = (hence 
is bounded), gives a bound for dY^d"^^^ A^/"^ , hence for dpd'l^'^^A. 

a: trivial 

ly; Apply i9™i9" on (p6|), then the right-hand side of ( [2^ ) is bounded by assumption, 
from which the claim follows immediately. 
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N: Again applying dj^d^ on the lapse equation (^) yields an equation of the form ( pl|) 
for df^dJ^N" plus some already bounded term (by assumption and the already proven 
bound for d^dJ^A'). The boundedness of d^^d^N' follows then from Gronwall's 
estimate after integrating the equation for dJ^d^N" along 5^ starting at a point 



with dY'd^N' = 0. 



H: Differentiation of equation (27) gives an expression of the form (dt — i'dx)d^^dJ^H' 



Bid^^dJ^H' + i?2, with Bi,B2 bounded, where the assumptions and the already 
obtained bound for dY^d^'^^v have been used. Thus applying Gronwall's inequality 
yields a bound for d^d'^'^^H. 

K: Apply 5™ 5" on equation (2C), then the right-hand side is bounded by the previous 
estimates, hence bounding d'^d'^'^^K. 



Proof of lemma W. 



A: Applying 5™ 5" on (|2^) yields immediately a bound for d^dl^A, since the right-hand 
side is already bounded. 

a: The same argument with 3™ instead of 5™ 5" applied to equation ( p4[ ) works in this 
case to bound d'^d'^a. 

H: And again, on equation (p^), yields an estimate for d^^d^H. 



K: Apply d^d^ on equation (|25D, then the only terms on the right-hand side, not 
already known to be bounded are d^^d^N" and d^^d'^A". But inserting equations 
(P), respectively (||) for N" and A", one sees easily, that dY'd^N" and d^d^A" 
indeed are bounded by the right-hand sides of their equations, hence d^^^d'^K is 
bounded. 

N: Unfortunately there is no explicit equation for A^. So we have to apply the deriva- 
tive operator 9™"*"^ 9^ to equation ([2l|), therefore producing already estimated terms 
involving dY~^^d^ applied on the quantities treated above, as well as on $ (bounded 
by property (^5[)), but also a term involving dY~^^d^~^^A on the right-hand side. If 
this term turns out to be bounded, then the equation for d^~^^d^N" has the form 

d^+^d^N" = B- A-^A'd^+^d^N' 

+ A^dr^^d^Nl^^iH - Kf + K^ + 47r(p + tr 5)) - A^ 

with \B\ bounded. On 5\ Sf+^OJA^ attains its maximum, from which we can infer 
ar+^a^iV < I (1 - B/A^) U{H - Kf + K^ + 47r(p + tr S)) | 

I J max 

hence it is bounded from above. Similarily for the minimum: 

9r+^9^Af > I (1 - B/A"^) (i(i7 - Kf + K^ + 47r(p + tr 5)) | 

V J min 

which is bounded from below and we get \dY'^^d'"N\ < C. 

It remains to show the boundedness of d]!^~^^ d"^^^ A = d'^d'^'^^A. Inserting equation 
(p^) for A, we see on the right-hand side beside a bounded term the quantities 



d"^dJ^A" and d^dJ^'^^u'. Inserting equation (19) for A" and equation (^q) for u', we 



finally obtain bounds for d^^d^A" and 5™5""^^z/', hence for d^d'^'^^A, as desired. 
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V. The final estimate is straightforward: Apply the operator 9!^ to equation 
(p^), which yields immdediately a bound for d^'^^dJ^i' since the right-hand side of 
the equation is already bounded by the arguments above. 



3.4 Higher order estimates 

Here we prove the matter regularity condition (^) for Einstein- Vlasov, Einstein-Maxwell 
and the Einstein- Vlasov-Maxwell system, achieving the goal of this work. 

3.4.1 Collisionless matter 

The 3+1-split of the Vlasov equation (|9|) with respect to the symmetries reads 
= dtf+(NA-^^-,^)dJ 



(36) + {-A-'N'v^ + NKv' + iV^-^^' 1^!)^)!) A / 



N{A-^A'^^-iiH-K))v^^f , 



where := \/l + 5ijV^v^ and = 9^.2/ = d^sf = v^^^f — v'^-^f (by symmetry) has 
been used. 



The energy momentum tensor (10) associated with the Vlasov equation leads to the matter 
quantities 

p = j fv^dv 
j = I fv^ dv 





Sab = / fVaVb/v dv 



Now we investigate the matter regularity property (35). Since / is constant along the 
characteristics of the Vlasov equation, given an initial particle distribution /o on the mass 
shell over S, we get the matter distribution for each time t as a function f{t,y,w) on 
the mass shell over the leaf St, by f{t,y,w) = fo{Y{t,y,w),V{t,y,w)), where {Y,V) 
denote the characteristic curve of the Vlasov equation through the point (0, y, w) € P(o,y), 
(0,2/) GS. 

Thus the matter quantities p, j, S are bounded, provided the support of / remains bounded. 
Define now 

Pf{t) := {suplvj \ v G supp f{t,y,w) y{Q,y,w)eP^o,y)} ' 

then the matter quantities are bounded by C(l + Pf{t))^, hence it is sufficient, to control 
Pf{t) on]ti,0]. 

Since the characteristic curves (Y, V) are integral curves of X and all coefficients in the 



components of X are already bounded by proposition 3.4, the characteristics itself are 
bounded, Pf{t) < C as desired. 

Now we need to iterate this procedure. Assume, that all didl.T and d^d^.^^ for k+l = m+n 
are bounded. To bound the derivatives of the matter quantities of order m + n + 1, 
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differentiate the Vlasov equation m + n + 1 times (with respect to t and x only). This 
yields linear equations for d^dl^f, for all non-negative integers k = m -\- n + 1, of the 
form 



X{d^dU) = B , 

where B vanishes outside the support of /, which bounds the support of B. Moreover, B 
involves derivatives of order m + n + 1 of the quantities J- and A', N'. For derivatives 
of order at most m in t, we see the boundedness of |-B|,_by applying similar arguments 
as used in the proof of lemma 3^. Thus in this case d^d^f is bounded and Pq^qt j, since 



they have the same characteristics, hence the derivatives of the matter quantities of order 
m + n + 1, involving derivatives of order at most m in t are bounded. Moreover, due to 
the simple dependence of the characteristics on the frame variables v, all derivatives of / 
with respect to v are bounded, too. 

In view of this fact, we are able to bound all derivatives of order m + n + 1 of the matter 
quantities, by considering the Vlasov equation as an equation for 9™"^^ 5"^^/, for which 
the right-hand side is known to be bounded. 

Therefore the first part of property ( |35| ) holds. For the second part we must show, that 
the spacelike derivatives of T can be redistributed to timelike derivatives of But this 
has already been done in the proof of the first part of (^) , and we get the 

3.8. Theorem. 

Let {M,g,f) be a surface symmetric solution of the Einstein- Vlasov system, which 
possesses a symmetric Cauchy surface S with strictly negative mean curvature H < 
H < and H = H somewhere on S. 

Then all of the past ofT, admits a PMC foliation {St}, where t takes all values in 
the interval ] — oo,0] and H takes all values in ] — oo,H]. 



3.4.2 Maxwell field 

Let F denote the electromagnetic field. The symmetry simplifies F: Due to the symmetries 
F can be written relative to the orthonormal coframe {V^} introduced in section 3A as 

F = -e{t, x) V° A + b{t, x) A , 

since all other components are forced to vanish and F and A V^, A %^ remain 
invariant under the action of the symmetry group. 

To obtain an explicit form of the Maxwell equations F and *F are represented in the given 
coordinates by 
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The Maxwell equations (^) in vacuum are given by setting J" = 0. We get for the magnetic 
field 



^23=0 
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6'=0 
6 = 



Since does not depend on {t, x), we find da{{Aa)'^ 6) = for a = 0, 1, hence 

b = C{Aa)-^ 
b = Ce^ . 



For the electric field we find 


*F23 = 



(K-^ey= 



(K-^e) 







(Note, that *-Foi,2 vanishes, since nb contains no e^). Again, since does not depend on 
(t, x), we get da{{Aa)'^ e) = 0, for a = 0, 1, hence 

e = C{Aa)-^ 

e = CNA-^a-^ . 

Since N,A^^,a^^ are already bounded, the same is true for the electromagnetic fields e 
and b. 

The energy momentum tensor (^) takes the simple form 

/I 0\ 



0-100 
10 
\0 1/ 



and we get the matter quantities as 



P 

j = 



Svr Stt ^ ' 



The property (|35| ) posed in section |3^ is obviously fulfilled, hence we get the 
3.9. Theorem. 

Let {M,g,F) be a surface symmetric solution of the Einstein-Maxwell system with 
plane or hyperbolic symmetry, which possesses a symmetric Cauchy surface S with 
strictly negative mean curvature H < H < and H = H somewhere on S. 
Then all of the past ofT, admits a PMC foliation {St}, where t takes all values in 
the interval ] — oo,0] and H takes all values in ] — qo,H]. 

Note that the restriction to plane and hyperbolic symmetry is necessary, since the electro- 
magnetic energy momentum tensor does not obey the non-negative pressures condition, 
and therefore fails to satisfy the assumptions of proposition 3.4 in the case of spherical 
symmetry. 
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3.4.3 Charged particles 

Now consider the Vlasov-Maxwell system in {M,g). We obtain the coupled equations by 
modifying the uncharged Vlasov equation (56|) by adding the term 



e + A-^A' - N-'u{l + A-^A%) ^ 
on the right-hand side. The matter current J according to the energy-momentum tensor 



e = T + E 



p= 9(^60/60) 
j=-G(%o/ei) 



is of the form 

J = p^eo + fei = N-'pdo + {A-'j - N-'up)di . 

The homogeneous Maxwell equations remain unchanged, which yields as before 

b = C{Aa)-^ 
b = Ce^ . 

Since *JaP'y = \/\9\^pLaP'yJ^ ■, ^ud again for a = 0, 1, da{n~^e) = \e^\da{{Aa)'^ e) the 
inhomogeneous equations read 



(K-ie)-=-47r^^ji 



e'=-2A-^A'e + AnAp 

'e = - 2{Aa)~'^ {Aa) e + Att{Avp - Nj) 



To get a bound for e, note that the inhomogeneous Maxwell equations are of the form 
dae = (fa e + ipa ■ {'') , a = 0, 1 



with ipa,il^a bounded by proposition Integrating the second equation gives 

<C f e +C , 



1 ./ol J' 



since J^i p, < C by integrating equation (19) as shown in the proof of proposition 

3.4. Using then the inequality 



{t)< 



(0) + 



and applying Gronwall's inequality on e| we obtain 
< C . 

Now using the first equation, we get analogously 

\e{t,y) - e{t,x)\ 



ry 










J 




[ Vie+ [ Vi-© 


< c 




J X 











+ c<c , 



where |/c.i e| < C has been used. 
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Thus we have 



\e{t,y)-eit,x)\ < C 

=^ \e\<C 



e 

51 



< C 



as desired. 



Now property (|35[) is obtained by bounding appropriate derivatives of e: AppUcation of 
on the equation for e' shows the first part of property (|35| ) and apphcation of 
9™ 9" on the equation for e establishes the second part. This proves: 

3.10. Theorem. 

Let {M,g,f,F) be a surface symmetric solution of the Einstein-Vlasov-Maxwell sys- 
tem with plane or hyperbolic symmetry, which possesses a symmetric Cauchy surface 
S with strictly negative mean curvature H < H < and H = H somewhere on S. 
Then all of the past ofT, admits a PMC foliation {St}, where t takes all values in 
the interval ] — 00, 0] and H takes all values in ] — 00,^]. 

3.5 Improving the results 

The purpose of this section is to get rid of the requirement of strictly negative mean 
curvature on the Cauchy surfaces, which seems to be a rather technical restriction. Let us 
instead assume S to be a symmetric Cauchy surface in the surface symmetric spacetime 
{M,g) with matter obeying the dominant and strong energy condition as well as the 
non-negative pressures condition, such as the surface symmetric Einstein- Vlasov system 
(M, (7, /). Note, that we assume the non- negative pressures condition not only for the 
spherically symmetric case. Indeed we will see, that this energy condition is necessary 
for the arguments given in this section, thus we cannot apply them to matter involving 
electromagnetic fields. 



As pointed out in 3.2.1 there is an important difference between the possible types of 
surface symmetry and it turns out to be a good idea, to analyse the spherically symmetric 
case separately. 

3.5.1 Spherically symmetric spacetimes 

There are some results already obtained elsewhere, so we can exclude some special cases 
from our analysis. Namely for the Einstein- Vlasov system (as well as for the massless 



scalar field) it has been proven in |R4] and |BR|, that given an arbitrary symmetric 
constant mean curvature Cauchy surface, there exist a global CMC foliation with the 
mean curvature taking all real values. So we may assume, that the mean curvature H on 



S is not constant, and we immediately get from theorem 2.2 the existence of a local in 
time PMC foliation ]ti,t2[y<T, of a neighbourhood of S in M. 

There are some a priori estimates: 

• It is shown in |0, that r and are bounded. Theorem 2.1 in [Q then shows, that 
all timelike curves have finite length. 



Inserting the bound for r on the right-hand side of the estimate (pj4|), we get from 
(^) bounds for r^^ and M on any finite time interval. 
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Section III in [BR] shows, that in hght of the fact, that all timelike curves have finite 



lengths, the volumes of two arbitrary Cauchy surfaces Si and S2 are related by 

VoKSa) < Vol($]i) ^1 + Csup\H\^ 

This implies, that the volumes of all Cauchy surfaces 5" in M are bounded above by 
the volume of S and, since H is bounded on each finite time interval, interchanging 
the roles of S and S shows, that the volumes are bounded from below, too, on each 
finite time interval. 



If we denote the volume form of a PMC leaf St by and the volume of St by V{t) we 
have = and V{t) = J ^ dx (I'd dip = A7ra~^ f^i r^. 

With the bounds of V, V^^,r, we get now o, < C, and again using the bounds for 
the radius function we get A, A^^ < C, thus the first fundamental form of the leaves is 
bounded from above and below on each finite time interval. 

Moreover, as shown in section |3.3| , the estimate ( |3^ bounds |yl~'^A'| and \K\, so that the 
second fundamental form is also bounded (from above) on each finite time interval as well 
as l^'l = \A^\ \A-^A'\. 

Consider now the lapse equation (|2l|), written in the form 

{AN'Y = A^n(UH - Kf + + Mp + S)) -A^ , 



where the term in brackets is non-negative. Therefore, setting this term to zero, we get 
the estimate {AN')' > -C, hence, for arbitrary p,q € St: {AN'){p) - {AN'){q) > -C . 
Now choosing q to be a critical point of N , yields N'(p) > — C, and similarily choosing p 
as a critical point of N gives N'(q) < C. From this we see, that |A^'| is bounded, since p, 
respectively q is arbitrary on St- 



The difference to the situation considered in section 3.3 is, that we are no longer confined 
to a mean curvature, having everywhere fixed sign. The difficulty with the estimates done 
there is, that they provide no information about the behaviour of the lapse function, when 
H becomes zero. The estimates done here so far, do not rely on this fact and the uniform 
bound of A'^' shows, that either lapse remains finite or diverges uniformly to infinity. 

In order to prove global existence of a symmetric PMC foliation assume, that ]ti,t2[ is 
the maximal time interval of existence. Without loss of generality let us consider only a 
possible extension towards the past, where H decreases with decreasing t. Thus, we are 



looking for regular symmetric initial data for t = ti in the sense of definition 3.2. 
Then there are two cases: First, H(t) is everywhere positive or everywhere negative near ti , 
the arguments in section |3.3| apply (the fixed sign of H is enough to perform the estimates, 
whether H is positive or negative), extending the foliation and we get a contradiction to 
the maximality of the time interval, hence ti = —00. 

Second, limt^ti H{t) possesses no unique sign near ti. Then we have to prove proposition 



3.4 and the two subsequent lemmas 3.6 and 3.7. The discussion at the beginning of this 



section has shown, that some of these quantities are already bounded. The crucial step is, 
to find a bound for lapse and its derivatives. 



The idea is to reparametrize the foliation as has been done for a CMC foliation in |BR]. 



The effect of the reparametrization on lapse and shift was outlined in section |2.1| . Let us 
introduce a function r by 

T := t + / N{u,X())du , 
Jti 
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where J^* means lims^^j J^. It is well defined, since / N{'y{u)) du along the integral curves 
of the normals of the leaves measures the length of 7, which we already know to be finite, 
and |A^'| < C ensures the integrability of N{t,xo) < N{'y{t)) + C. This construction 
works in either case, whether is bounded or diverges uniformly to infinity. The function 
t I — > T{t) is monotone since ^ = 1+N{t, xq) and turns out to be an orientation preserving 
diffeomorphism ]ti,i2[ — * ]ti = ti,T2[, thus r can be used as a new time function for 
the foliation, which squeezes the time differences of neighbouring leaves compared with 
the normal vector by adding to t the normal component of the length of the piece of the 
xo = const line connecting {ti,XQ) with {t,XQ). This produces a new lapse function and a 
new shift vector: 

N 



\dt J l + N{t,xo) 

"'ydtj ''~l + N{t,xo) 

Conversely, t stretches time differences of adjacent leaves, such that we get the inverse 
transformation by i = t — J^^ ^{u, xq) du, which subtracts from r the former added 
'length' of the xq = const line now expressed in terms of the new time coordinate. By 
the inverse transformation, we see that ^ = 1 — N{t,xo), and therefore get the relation 
1 — N{t, xq) = (1 + N{t, a^o))""^ in points with spatial coordinate xq. 

The benefit of this reparametrization is, that due to |A^'[ < C, N and |A^'| remain bounded, 
and we could try to analyse the field equations according to the reparametrized foliation. 
Inspection of these equations shows, that we get the field equations in the new coordi- 
nates from equations ([l9|)-(|27|) by replacing everywhere dt by 9,-, N hy N, 1/ hy u (since 
reformulating of the constraint and evolution equations preserves their form) while the 
lapse equation (|2l| ) and the PMC condition (^) depend on the parametrization, and 
hence have to be modified: The subtraction of on the right-hand side stems from the 
term A^Nn{H), where the PMC condition sets Nn{H) equal to one. Replacing N hy N 
yields {Nn{H)){l + N{t,xo))-^ = (1 + N{t,xo))-^ = 1 - N{t,xo). The PMC condition 



(22), which holds in the old parametrization, can be formulated in the new coordinates by 
expressing u in terms of z>, which produces an additional summand 1 — N{t,xq) on the 
right-hand side, besides the replacement of v by D. 

Now we start to analyse the new field equations, trying to retain the same line of thought 

A'\,a,a 



Hi K 



as well 
n 



of section 3.3: The analysis already done here bounds A,A~ , 

as, of course, N and \N'\. The arguments carried out in section bound now and 
\H'\. Finally we need a bound for N~^. Since it is not clear, in the case where A'^ diverges 
near ti, whether N grows monotonically or not, we divide the time interval ]ii,t2[ into 
two subsets. First consider all points, where N{t, xq) > M, with some suitable chosen real 
number M, which will be specified later. Then we get the estimate 



1 - N{t,x) 



1 



N{t,x) 



l + N{t,xo) 



1 + N{t,xo) - N{t,x) 



l + N{t,xo) 



< 



l + C 

1 + M 



where C is an upper bound for \N{t,xo) — N{t,x)\, whose existence is guaranteed by the 
bounds for lA^'l and the volume of the leaves of the foliation. Choose now M = 2C and 



we get |1 - iV| < 



i+c 

1+2C 



< 1, thus A'^ is bounded away from zero for all points, where 



^itjXo) > M. Consider now the points, where N(t,xo) < M. Then we have 



1-N{t,xo: 



N{t,xo) 
l + N{t,xo) 



l + N{t,xo)-N{t,xo) 
l + N{t,xo) 



> 



1 



1 + M 
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Now we can apply the arguments in p5[| to the modified lapse equation, as has been done 
in section (whereupon the factor 1 — N{t,xq) causes no trouble, since it is bounded 
from above and below) and we get a bound for A^~^ for all points, where N{t,XQ) < M, 
and we are done. 



Therefore we have estimated all quantities appearing in proposition 3.4. lemma 3.6 is also 
true for the new field equations without the need to modify the proof, while the proof of 
lemma |3.7| has to be modified, since the argument bounding \d^~^^d!^N\ does not carry 
over. The lapse equation and the equation for its time derivatives are 

{AN')' = A^n(^^{H - K)^ + K^ + Mp + tr 5)) -A^(l- N{t, xq) ] 
{Ad^N')' = A^dTN{\{H - K)^ + K^ + 4tt{p + tr 5) ' 
-A'(l-drN{r,xo))+B , 
where B denotes an already bounded quantity at the corresponding stage in the inductive 



argument in lemma 3.7. In this situation lemma 1 in ||BR| applies literally, bounding the 



time derivatives of the lapse function, therefore completing the proof of lemma 3.7. 
Putting all arguments together we are able to extend the foliation beyond ti, contradicting 
the assumed maximality of the interval of existence, thus arriving at the following 

3.11. Theorem. 

Let {M,g,f) be a surface symmetric solution of the Einstein- Vlasov system with 
spherical symmetry. 

Then we can foliate the whole spacetime by a PMC foliation, where the time function 
takes on all real values and the mean curvature of the leaves tends uniformly to ±oo 
for t zboo, thus producing crushing singularities. 

This theorem provides barrier surfaces (see [0), establishing now the existence of CMC 
hypersurfaces for each value of the mean curvature, therefore proving the closed universe 
recollapse conjecture in this case: 

3.12. Corollary. 

In the situation of theorem 3.11 , the spacetime possesses a global CMC foliation and 
the mean curvature takes on all real values. In particular, the spacetime admits a 
maximal Cauchy surface. 



3.5.2 Spacetimes with plane and hyperbolic symmetry 

We will see, that most of the arguments performed in the spherically symmetric case carry 



over to the past domain of dependence D~{Ti) in the expanding models (compare 3.2.1 for 
a precise definition of this terminology). Due to this fact we perform the following analysis 
only on the half-open time interval ]ti,0] and assume further, {M,g) to be non-flat in the 
plane symmetric case and the mass function (|3T| ) to be positive on S in the hyperbolic 
case. 

Again we can exclude some cases, already investigated (although this does not matter. 



since the arguments here apply to these cases). It has been proven in [114] and pR| , that 
given a symmetric Cauchy surface with negative constant mean curvature Hq (remember 
the restrictions on the mean curvature in the expanding models, compare |3.2.1[ ), there 
exists a global CMC foliation with the mean curvature taking all values in the interval 
] —oo,Hq]. So we can assume without loss of generality, that the mean curvature H on T, 
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is not constant and not everywhere positive. Again, this assumption ensures the existence 
of a local in time PMC foliation x^] of a neighbourhood of S in M . 

Now we find similar a priori bounds 

• In Z)^(S) r is bounded, since dr is future pointing timelike everywhere and S is 
compact. is bounded on S by assumption in the hyperbolic case and in the 

plane symmetric case this follows from lemma 2.4 in |R4]. By the mass flux equation 



(p3|) together with the non-negative pressures condition, M grows along past pointing 
timelike curves, thus is bounded in all of Z)~(S). Theorem 2.1 in Q then shows, 
that all timelike curves in have finite length. 



• Applying the arguments in the corresponding place of section |3.5.1| we get bounds 
for and M for any finite time interval of the form ]ti,0]. 

• Again, the corresponding argument performed in the spherically symmetric case 
holds, and we get upper and lower bounds for the volumes of arbitrary Cauchy 
surfaces in M. 

Finally, the volume of the leaves St is given by V{t) = Ca^^ f^i r^, establishing bounds 



for the first fundamental form and its inverse (in D (S)) as has been shown in 3.5.1 
therefore all of the remaining arguments performed there apply and we get the 

3.13. Theorem. 

Let {M,g,f) be a surface symmetric solution of the Einstein- Vlasov system with 
plane or hyperbolic symmetry and S a symmetric Cauchy surface in M . 
If {M,g) is non-flat in the plane symmetric case and the mass function is positive 
on S in the hyperbolic case, then we can foliate all of the past of S by PMC hyper- 
surfaces, where the time function takes on all values in the interval ] — oo, 0] and the 
mean curvature of the leaves tends uniformly to — oo for t — > —oo. 

Using again the PMC leaves as barrier surfaces we get the 

3.14. Corollary. 



In the situation of theorem 3. IS D (S) possesses a CMC Cauchy surface for each 



value of the mean curvature in ] —oo, min^ II[. 



3.5.3 Comparing the results 



To close this chapter, there are some remarks concerning the theorems 3.11 and 3.13| and 
the relation to the work done in the preceding sections. 



At first, it is obvious, that theorem 3.11 generalizes theorem 3.8 in the spherically symmet- 



ric case. In the plane and hyperbolic case theorem 3.13 does not generalize the theorems 



|, |3.9| and |3.1C| . The latter ones only establish a global PMC foliation unless the mean 
curvature of the leaves becomes zero somewhere (and again the leaves provide barrier 
surfaces for CMC Cauchy surfaces, establishing corollary |3.14 in the situation of those 
theorems). Theorem 3.13| is not restricted to this condition, but it needs the extra non- 
negative pressures condition, which excludes electromagnetic fields. 

Apart from the non-negative pressures condition the assumption of positive mass on S is 
a non-trivial constraint in the hyperbolic case, while automatically fulfilled in the non-flat 



plane symmetric case by lemma 2.4 in [R.4]. As shown above, the positivity of mass is 
needed to obtain a bound for the length of timelike curves in M, by applying the arguments 
, thus necessary to the construction done here. 



m 
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4 Conclusion and outlook 



First I list the main results achieved in this work 

• For the spacetimes considered so far the existence of a global PMC foliation has been 
shown, where several matter models have been taken into account. 

• For the model, which is not expanding or contracting everywhere, the closed uni- 
verse recollapse conjecture has been proved. In particular the foliation covers the 
whole Cauchy development of the initial hypersurface, with a crushing singularity 
both in the distant past and future. Moreover the spacetime possesses a maximal 
hypersurface. 

• In the expanding models the foliation covers at least the whole past of the initial 
hypersurface towards a crushing singularity. 

The choice of matter turned out to be important only in so far, as some energy conditions 
are satisfied and the matter fields do not become singular in a regular geometric back- 
ground. 

The necessary energy conditions are the dominant and strong energy conditions. Further- 
more to obtain stronger results, which do not rely on a fixed sign condition for the mean 
curvature, the non-negative pressures condition is required, to ensure that the Hawking 
mass does not tend to zero, contributing enough attraction, that the lifetime of timelike 
curves become finite. 

The rather strong results about locally spatially homogeneous spacetimes in [ |R3|| likewise 
rely on this condition, which appears in the slightly relaxed form, that only the sum of 
the pressures has to be non-negative, due to the high degree of symmetry in those models. 
This relaxed condition even permits electromagnetic fields, a type of matter, which does 
not satisfy the non-negative pressures condition, thus not leading to the stronger results 
in the context of spacetimes with less symmetry. 

The close analogy of the proofs for the global existence of CMC and PMC foliations indi- 
cates, that all results obtained for CMC foliations may also be proved for PMC foliations. 
If this conjecture turns out to be true (the second part of this paper will give another pos- 
itive example for this conjecture), then we are no longer concerned with the topological 
restrictions imposed on the existence of CMC hypersurfaces, having a much more fiexible 
tool at hand. Moreover, the global results on PMC foliations can be used to provide bar- 
rier surfaces, which guarantee the existence of a CMC foliation, where the mean curvature 



ranges between these barriers (compare the corollaries 3.12 and 3.14]) . 

The results demonstrate, that in the cases considered in this work a satisfying answer has 
been given to the global existence problem for PMC foliations. Then the related question 
arises, whether the foliation covers the whole spacetime. The present work gives only a 
partial answer to this question. Denoting the initial Cauchy surface by S, we saw, that we 
have covered the whole past Z)~(S) in the expanding spacetimes and the whole Cauchy 
development -D(S) in the recollapsing models. These positive results have two limits 

1. In the expanding models there remains an open question about the future of T, in 
M. Either the mean curvature becomes zero somewhere or there is no control on the 
radius function towards the expanding direction. In each case the present techniques 
do not apply. In addition there are topological obstructions in the expanding direc- 
tion, preventing the mean curvature to become positive, a foliation of the future of 



30 



T, must stop before this happens. But this does not imply, that the whole of D~^{Y,) 
can be covered by such a foliation, since the leaves of the foliation may become null 
or non-compact, where the general notion of singularity avoidance comes up. 
Namely, compared with the strong results obtained in pl3[ in the locally spatially 
homogeneous models, we require more information about the asymptotic behaviour 
of the foliation to obtain results on geodesic completeness. 

2. In the contracting direction and in the recollapsing models things look different. In 
the contracting direction the theorems ensure the existence of crushing singulari- 
ties, and all of the Cauchy development of S is covered by the foliation. This is 
a consequence of Hawking's famous incompleteness theorem for globally hyperbolic 
spacetimes, satisfying the timelike convergence condition, where the maximal time 
of existence is estimated by 3/\H\ and H — > — oo. But there is no information about 
the boundary of D{Y,), which might be either a curvature singularity or merely a 
Cauchy horizon. 

The work in pi3| ] indicates roughly, what remains to do. Relating the particle cur- 
rent density to the energy density and investigating the asymptotic behaviour of the 
former may give rise to unbounded curvature and produce a curvature singularity. 

For the spacetimes considered here, the work p^e| of Rein is important. He succeeded in 
finding satisfying answers to the above questions, by using a time function intimately tied 
to the symmetry. Unfortunately it is not clear how to generalize his approach to other 
spacetimes. Although the estimates done here also exploit the high symmetry of the mod- 
els, the construction itself does not depend on it, motivating this work. A generalization 
to spacetimes with two commuting local Killing vector fields will be the content of the 
second part. 



31 



References 



[A] H. Amann, Gewdhnliche Differentialgleichungen, Walter de Gruyter, 1995. 

[An] Hakan Andreasson, Global foliations of matter spacetimes with Gowdy symmetry, 
Commun. Math. Phys. 206 (1999) 337-365 

[B] G.A. Burnett, Incompleteness theorems for the spherically symmetric spacetimes. 
Physical Review D 43, No. 4 (1991), 1143-1149 

[BR] G.A. Burnett, A.D. Rendall, Existence of maximal hypersurfaces in some spherically 
symmetric spacetimes, Classical and Quantum Gravity 13 (1996), 111-124 

[Ba] R. Bartnik, Remarks on cosmological spacetimes and constant mean curvature sur- 
faces, Commun. Math. Phys. 117 (1988), 615-624 

[FR] H. Friedrich, A.D. Rendall, The Cauchy problem for the Einstein Equations in B.G. 
Schmidt (ed), Einstein's Field Equations and Their Physical Implications, Lecture 
Notes in Physics 540, (2000) 

[G] C. Gerhardt, H-Surfaces in Lorentzian Manifolds, Commun. Math. Phys. 89 (1983), 
523-553 

[GM] J. Guven, N. 6 Murchadha, The constraints in spherically symmetric classical gen- 
eral relativity I, Physical Review D 52 (1995), 758-775 

[HI] O. Henkel, Prescribed Mean Curvature foliations in cosmological spacetimes with 
matter, Dissertation, Universitat Potsdam, 2001 

[H2] O. Henkel, Local Prescribed Mean Curvature foliations in cosmological spacetimes, 
gr-qc/0108003 

[HE] S.W. Hawking, G.F.R. Ellis, The large scale structure of spacetime, Cambridge Uni- 
versity Press 1973 

[IR] J. Isenberg, A.D. Rendall, Cosmological spacetimes not covered by a constant mean 
curvature slicing. Class. Quantum Grav. 15, (1998), 3679-3688 

[Re] G. Rein, Cosmological solutions of the Vlasov-Einstein system with spherical, plane, 
and hyperboHc symmetry. Math. Proc. Camb. Phil. Soc. 119 (1996), 739-762 

[Rl] A.D. Rendall, Constant mean curvature foliations in cosmological spacetimes, Helv. 
Phys. Acta 69 (1996), 490-500 

[R2] A.D. Rendall, An introduction to the Einstein- Vlasov system, Banach Centre Pub- 
lications 41 (1997), 35-68 

[R3] A.D. Rendall, Global properties of locally spatially homogeneous cosmological mod- 
els with matter. Math. Proc. Camb. Phil. Soc. 118 (1995), 511-526 

[R4] A.D. Rendall, Crushing singularities in spacetimes with spherical, plane and hyper- 
bolic symmetry. Classical and Quantum Gravity 12 (1995), 1517-1533 

[R5] A.D. Rendall, Existence and non-existence results for global constant mean curvature 
foliations. Nonlinear Analysis, Theory, Methods and Applications 30 (1997), 3589- 
3598 



32 



A.D. Rendall, Existence of constant mean curvature foliations in spacetimes with 
two-dimensional local symmetry, Commun. Math. Phys. 159 (1997), 145-164 

R. Schoen, S.-T. Yau, Existence of incompressible minimal surfaces and the topol- 
ogy of three dimensional manifolds with non-negative scalar curvature. Annals of 
Mathematics 110 (1979), 127-142, 

R.M. Wald, General Relativity, The University of Chicago Press 1984 



33 



